Scale of Homogeneity of the Universe from WMAP 
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We review the physics of the Grishchuck-Zel'dovich effect which describes the impact of large ampli- 
tude, super-horizon gravitational field fluctuations on the Cosmic Microwave Background anisotropy 
power spectrum. Using the latest determination of the spectrum by WMAP, we infer a lower limit 
on the present length-scale of such fluctuations of 3.9 x 10 3 times the cosmological particle horizon 
(at the 95% confidence level). 
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A definitive map of large scale fluctuations in the Cosmic 
Microwave Background (CMB) is currently being ren- 
dered by the WMAP satellite [1]. The preliminary re- 
sults have been used to constrain a host of cosmological 
parameters within the context of a Big Bang universe 
with primordial adiabatic perturbations [2,3]. A candi- 
date model, which consists of a cocktail of cold dark, 
matter, baryons, neutrinos, radiation and a cosmological 
constant (the ACDM model), seems to fit the data well, 
albeit with a larger goodness of fit than would be de- 
sired. The paradigm of an essentially homogeneous uni- 
verse with mild deviations is systematically passing many 
of the observational tests with which it is confronted. 

The homogeneity of space-time is one of the funda- 
mental assumptions of our current working model of the 
universe. One possibility is that our observable universe 
is in a particularly quiet patch of an inhomogeneous space 
time which has undergone a period of inflation [4-6]. 
Conceivably there may be indirect evidence of the pri- 
mordial inhomogeneity within which our homogeneous 
patch resides. Grishchuck and Zel'dovich (GZ hereafter) 
have pointed out that such large scale inhomogeneity will 
have an impact on the amplitude of the quadrupole of 
the CMB [7]. In this report we combine their formal- 
ism with a goodness of fit analysis of the WMAP data to 
find a lower bound on the scale of inhomogeneity of the 
universe. 

On large scales, the Sachs- Wolfe (SW) effect [8] re- 
lates metric perturbations with the CMB temperature 
anisotropics. It is given, along a direction n in the sky, 
in terms of the perturbations in the gravitational poten- 
tial, $, by 
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where r\ is the conformal time defined today as r\§. The 
conformal time today is the comoving distance to the 
horizon. The second term is the Integrated Sachs- Wolfe 
(ISW) effect and accounts for the contributions to the 
fluctuations due to the possible variations of the gravita- 
tional potential along the line of sight whereas the first 
term is due to perturbations in the potential present at 
the last scattering (LS) surface. 



By expanding the temperature fluctuations in spherical 
harmonics on the sky, it is straightforward to show that, 
in a flat Universe, 
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where we have used ^(x, rf) = J $(k, 77)e~ lxk G? 3 fc/(27r) 3 
and k is the comoving wave-number related to x by 
k = 2ir/x. For the suppression factor <?, defined as 
<I>(k, 77) = g(r))$(k, t]ls), we use the precise parametriza- 
tion of [10]. Note that g(r)Ls) = 1 an d g varies with time 
depending on the cosmological model considered. P^ s is 
the gravitational field power spectrum at the last scatter- 
ing surface. With a flat spectrum P^ s (k) = Ak~ 3 , scales 
much larger than the horizon make little contribution to 
the Cis: If we consider only the SW term in Eq. (2), 
the wave-number contributions such that k <C 1 /r?o scale 
with k 21 . But a possible sharp and localized increase in 
the power spectrum at a certain super-horizon scale kaz 
of the type P$ s (k) = A5(k — koz) would leave a stronger 
imprint in the C/s at low multipoles as firstly shown 
by Grishchuck and Zel'dovich in [7]. This effect is more 
prominent for the monopole, dipole and quadrupole com- 
ponents. The monopole cannot be distinguished from 
a small shift in the measured temperature of the CMB 
because the density wave is super-horizon sized and we 
are averaging over our horizon. The dipole is contam- 
inated by the Doppler effect caused by the motion of 
our Local Group with respect to the CMB rest-frame. 
The relevant multipole to study the GZ effect is there- 
fore the quadrupole. For this multipole, contributions 
from super-horizon gravitational inhomogeneous modes 
of characteristic scale k to the GZ effect get suppressed 
by a factor of 0(k 4 ), as will be shown later. By not seeing 
an increase of power in the quadrupole, one can therefore 
place constrains on the homogeneity properties of super- 
horizon perturbations namely on the minimum allowed 
wavelength. On larger scales, due to the suppression fac- 
tor of 0(k ), the universe can be anisotropic and inhomo- 
geneous without being in disagreement with the available 
experimental constraints [9]. In fact, if we were to con- 
sider contributions from larger scales as well, (assuming, 
for example, a step function P£ s (k) — A9(koz — k)), 
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the total contribution to the quadrupole value would de- 
crease, doubling the value obtained for k GZ from the 
data. This behaviour is due to the way the normaliza- 
tion constant A is derived. As will be explained below, 
in this case, when fixing A, we would also be considering 
contributions from modes smaller than k GZ which would 
reduce the relative amplitude of the potential fluctua- 
tions coming solely from the k GZ mode, as compared to 
the standard Grishchuck and Zel'dovich calculation. 

We consider that the homogeneity of the gravitational 
field fluctuations is satisfied as long as the characteristic 
amplitudes of the fluctuations are below one (in the lin- 
ear regime boundary) . In the limiting case of amplitudes 
of order 1 , which we will investigate, the amplitude A of 
the power spectrum can be obtained by imposing that 
< |<J>(L GZ )| 2 >= 1, where < \§(L GZ )\ 2 > is the average 
variance of the gravitational field fluctuations in spheres 
of comoving radius L GZ /2 given by L GZ /2 = 2-n/k GZ . 
We point out that there is some confusion in the litera- 
ture in this respect [13,14]: we are fixing the amplitude of 
the gravitational field perturbation associated with the 
density perturbation, not the amplitude of the density 
perturbation. 

The normalized expression for the Grishchuck- 
Zel'dovich effect is the following 
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At I = 2, for kaz — 10~ 6 Mpc -1 and for the cosmologies 
considered here, the effect of the ISW corresponds to 
around 10 % of the pure SW contributions. Looking at 
the purely SW term of this expression, which dominates 
over the ISW term, we can use the small- argument limit 
of the spherical Bessel function J2{x) = x 2 /15, to obtain 
the following relation for the quadrupole amplitude 

^L(l = 2)^(k Gzm f = (2^) 2 (4) 

where AT/T = y^+IjQA^r), L = 2 m is the 
present size of the horizon diameter and L GZ is the 
present length-scale of the super-horizon homogeneous 
patch of characteristic wave-number kaz- This approxi- 
mation was obtained previously by others [9] without the 
2ir factor which was neglected. Having an upper limit on 
the value of the quadrupole, we can obtain a lower limit 
on the diameter size of the largest possible homogeneous 
scale of the universe. Similar expressions for the GZ ef- 
fect were calculated for the open cosmology case [13] for 
which the constraints on L GZ can be even more strin- 
gent. Given the present experimental consensus around 
the flatness of our universe, we place ourselves in a flat 
geometry. 

We proceed by using the recent WMAP temperature 
data to constrain the possible contribution of super- 



horizon gravitational field fluctuations (with an ampli- 
tude of order 1) to the quadrupole. For the analysis 
we consider the GZ contribution to the power spectrum 
in addition to the usual power law A-CDM spectrum. 
We compute the standard Ce spectrum by using CAMB 
code [11], and add to it the GZ spectrum computed an- 
alytically as described in Eq. (3). We consider a grid of 
models corresponding to inflationary adiabatic perturba- 
tions in a flat cosmology by varying the following cosmo- 
logical parameters around the best values found by the 
WMAP collaboration [2]: iJ , ^a, t, n and L GZ /L a . 
The normalization of the "power law" Ci spectrum was 
left free and was marginalized over. No contributions 
from neutrinos or gravitational waves were considered. 
As the effect of the super-horizon perturbation is ex- 
pected at very low I (mainly at the quadrupole), where 
the presence of the baryons is negligible, we fixed the 
baryon contribution at the best estimate from WMAP: 
il b h 2 = 0.023. We used the likelihood code provided by 
the WMAP team [3] and modified by [12] to compute 
the likelihood of each model of our grid. We then maxi- 
mized over all the parameters, except L GZ /L , in order 
to retrieve the likelihood as a function of the scale of 
the super-horizon perturbations. The result is plotted in 
Fig. 1. 

We obtain an upper limit on the region of homo- 
geneity (and isotropy) of L GZ > 3927Lo (at 95% CL). 
This correspond to L GZ ~ 6.9 x 10 7 Mpc for H n = 68 
Kms _1 Mpc _1 . As expected this value is comparable (ex- 
cept for the 2tt factor, as detailed above) to previous 
constraints relying on the COBE data [13-15], which also 
presented a low quadrupole value. The estimate from [16] 
gives Qcobe = 10±3±7 [iK, where the quoted errors are 
statistical and systematic (from the foreground removal) 
respectively. The WMAP results [17] give roughly the 
same upper limit of Qwmap = 12.3±3.1 [iK at 68% CL, 
where the systematics were assumed to be negligible. As 
stressed before, the contribution from the GZ effect to the 
quadrupole is proportional to k G z such that observing a 
quadrupole amplitude higher by a factor of 2 decreases 
the length-scale L GZ by only a factor of 2 1 / 4 . The small 
difference in the upper limit for the quadrupole between 
the WMAP and the COBE experiments induces an even 
smaller difference on the determination of the possible 
scale of the gravitational fluctuations considered. Never- 
theless our result is more robust because we include the 
present uncertainty on most of the relevant cosmological 
parameters. It is worth mentioning that - as could have 
been predicted - nearly no degeneracies were found be- 
tween the scale of the super-horizon fluctuations and all 
the other parameters explored in this analysis. This is 
due to the unique and strong localization of the GZ effect 
at the quadrupole component. 

There has been some debate surrounding the statisti- 
cal significance of the low quadrupole values obtained by 
the WMAP and the COBE DMR experiments as com- 
pared to the values predicted by the standard cosmo- 
logical model [2,17-21]. In this context, some work has 
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FIG. 1. Likelihood of the ratio Lgz/Lq. Lgz is the 
length-scale of super-horizon gravitational field fluctuations 
of amplitude of order 1 inferred by the WMAP temperature 
data. Lo is the horizon scale. All the other cosmological pa- 
rameters investigated have been maximized over. The lower 
horizontal line shows the 95% and the higher one, the 68% 
confidence level on Lgz /Lo. 

been devoted to finding a physical mechanism able to 
explain the possible discrepancy [2,19,18,22-29]. In this 
report, we take the quadrupole and octopole results from 
WMAP at face value and assume they are in accordance 
with the standard model within the cosmic variance er- 
rors, regardless of the present controversy. Of course, if 
any of these proposed mechanisms, modifying the stan- 
dard scenario, are indeed lowering the quadrupole value 
than our constraints on the homogeneity scale would be 
slightly affected. 

To conclude we point out that one can estimate the 
number of e- foldings, N e , which a patch of the universe 
has undergone during an inflationary period. Numer- 
ical studies of the onset of inflation in inhomogeneous 
universes indicate that the region of space time which 
inflates must start of with at most fluctuations of order 
one on the horizon scale [30,31]. Mapping this scale onto 
L GZ today we find that N e ~ ln[M in/ /(10" 15 GeV)] 
where M in f is the energy scale of inflation. For inflation 
occurring at M inf > 10 16 GeV one has N e > 72. It is 
interesting (but not surprising) that N e is similar to the 
number of e-folding required to get the correct amplitude 
and slope of fluctuations from the parametric amplifi- 
cation of quantum fluctuations during inflation [32,33]. 
One can also reinterpret this result as a constraint on 
the duration of an uninterrupted period of inflation and 
in doing so limit the effects of primeval, non-thermal cor- 
rections to cosmological perturbations [34] . 
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